Site determination and thermally assisted tunneling in homogenous nucleation 
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A combined low-temperature scanning tunneling microscopy and density functional theory study 
on the binding and diffusion of copper monomers, dimers, and trimers adsorbed on Cu(lll) is 
presented. Whereas atoms in trimers are found in fee sites only, monomers as well as atoms in 
dimers can occupy the stable fee as well as the met ast able hep site. In fact the dimer fcc-hcp 
configuration was found to be only 1.3 meV less favorable with respect to the fcc-fcc configuration. 
This enables a confined intra-cell dimer motion, which at temperatures below 5 K is dominated by 
thermally assisted tunneling. 

PACS numbers: 68.35.Fx, 68.37.Ef, 68.55.Ac, 66.35.+a 



In recent years important progress has been made in 
understanding epitaxial growth processes on the atomic 
scale, from both the experimental and the theoretical 
side. Early experiments employed field ion microscopy 
(FIM), whereas in recent years scanning tunneling mi- 
croscopy (STM) proved to be a perfect technique for this 
task. 

Interestingly, the first steps in the homoepitaxial 
growth on fee (111) noble- metal surfaces still belong to 
the unsolved problems. Compared to the more open 
metal surfaces, the behavior of single adatoms is more 
complicated on (111) surfaces, because adatoms can re- 
side in two different binding sites, the so-called fee and 
hep sites. Whereas the fee site continues the ABC layer 
stacking of the bulk the hep site breaks this sequence. 
As the energy difference between the two sites and the 
related diffusion barriers can be very small, the site oc- 
cupation preference is expected to depend critically on 
the size of the metal structures, i.e., it is different for 
monomers, dimers and trimers as was demonstrated ex- 
perimentally for Ir/Ir(lll) [ll0. 

The experimental determination of the fcc/hcp site 
preference of individual adatoms is a difficult task, and 
only two experimental studies exist. Wang and Ehrlich 
determined single adatoms of Ir/Ir(lll) to be located on 
hep sites 0,3, whereas in the case of Pt/Pt(lll) a clear 
fee site preference was found 5] . A general rule predicting 
the site preference was later developed for metals having 
a partly filled d-band |||, |7] that can be applied in both 
cases. Whereas no experimental work can be found in the 
literature for systems without a partly filled (i-shell, sev- 
eral theoretical studies exist. First-principles theoretical 
work by Stumpf and Scheffier |^] and Feibelman 9] pre- 
dicted an hep site preference for individual Al atoms on 
Al(lll). Other theoretical studies on copper, silver, and 
gold do not find any clear hep or fee site preference for a 
single adatom adsorbed on the fcc(lll) surface [Tol ITU . 
Clearly, very subtle contributions of both the electronic 
structure and the local relaxation have to be taken into 
account ^3 5 therefore no general rule is available to 
predict the preferred binding site. Although in homoepi- 



taxy Ir/Ir(lll) is the only known experimental example 
of an hep site preference so far, several similar systems 
may well exist. 

In this work we determined the adsorption site of Cu 
metal adatoms to be the fee site, by means of a unique 
technique of STM diffusion studies in combination with 
atomic manipulation, applicable also for systems that 
cannot be studied by FIM. Moreover a particular local- 
ized diffusion behavior of Cu dimers has been observed 
experimentally, with the properties as predicted earlier 
in a theoretical work by Bogicevic et al Despite 
the large mass of Cu atoms, this localized diffusion is 
dominated below 5 K by thermally assisted tunneling of 
Cu adatoms. Copper as a noble metal is characterized 
by an almost complete d-band occupation, making the 
existing, qualitative theoretical analysis of site prefer- 
ence ^ indeterminate. Instead the experimental re- 
sults are complemented by an extensive set of large-scale 
ab-initio density functional theory (DFT) calculations of 
the site preference on Cu(lll). 

Our experiments were performed with a low- 
temperature STM operated at 4 — 21 K. The sample was 
cleaned by sputtering and annealing cycles. Bias voltages 
refer to the sample voltage with respect to the tip. To 
measure possible tip effects we determined the hopping 
rate for different tunneling parameters and varied the 
tip-adatom interaction time by taking image series with 
fixed parameters but for different time intervals. Using 
atomic manipulation the monomers and dimers were po- 
sitioned far away from defect sites. Great care was taken 
to calibrate the temperature measurements precisely. 

Our calculations are within the generalized gradient 
approximation (GGA) using the plane-wave ultrasoft- 
pseudopotential code DACAPO 14]. We used slab su- 
percell geometries with a vacuum separation ^ 13 A, and 
with a five- (four-)layer slab thickness and a 3-by-3 (4- 
by-5) atom in-surface extension to map the monomer 
(dimer) adsorption and dynamics. The two lowermost 
layers were kept frozen in a bulk-lattice structure. The 
adsorbates and surface atoms were relaxed to 0.05 eV/A 
using an energy cutoff of 30 Ry, with dipole correction en- 
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abled. Our DFT calculations should determine the local 
relaxation, binding, and vibrations extremely accurately. 
However, ab-initio DFT can encounter difficulties in re- 
solving the energy variation between different adsorption 
sites if they are in the range of only a few meV 15]. Con- 
sequently, we have tested the robustness of our results 
and used a carefully controlled sampling of /c-points: a 
standard choice of 2 x 2 with a 0.1 eV Fermi smearing 
tested against 4x4 /c-point calculations with a 0.025 eV 
smearing. In addition, we have tested the monomer site- 
preference result against calculations using a fully relaxed 
four-layer 5-by-4 unit cell. 

The diffusion paths of individual adatoms were 
mapped by taking a series of consecutive STM images 
on Cu/Cu(lll). Five series of more than 100 images 
each were taken at different temperatures between 1 1 and 
12 K, allowing the determination of the diffusion barrier 
as Eb = (37 ± 5) meV 16j, in good agreement with other 
experiments XLilM ^^^^ theory 10, 19]. The map- 
ping of the relative movement of one adatom in between 
subsequent images clearly reveals the hexagonal lattice of 
the substrate surface, showing that only one adsorption 
site is observed. As pointed out by Wang and Ehrlich [j, 
the conventional technique of mapping diffusion paths is 
not applicable for a site determination if only one of the 
sites is occupied. On the other hand, lateral manipula- 
tion on the atomic scale [23, 0| can in general be used 
as a tool to access also metastable sites by artificial oc- 
cupation of non-equilibrium configurations ^22]. 

For this type of extension of the mapping technique, 
individual adatoms were repeatedly repositioned using 
lateral manipulation at a temperature of 7 K. At this tem- 
perature the monomers can be positioned in both the fee 
and the hep sites, making the site determination possible. 
If the temperature is increased to 10 K the individual 
adatoms regain their equilibrium position, which can now 
be identified as the fee site (see Fig.QJ. 

In order to have a point of origin for this extended 
site-mapping technique, the positions were determined 
with respect to a monatomic copper row built up atom 
by atom. The row ran parallel to one of the close-packed 
[110] directions. The lateral distances between the in- 
dividual adatoms and the row were \Ax\ = n • Y^3/8ao, 
with integer values of n for the fee positions and frac- 
tional values of n for the hep positions. In addition, this 
shows that the atoms in the row reside in fee positions 
in agreement with our theoretical calculations 23]. 

The diffusion experiments, in which no hep monomer 
occupation could be identified, revealed that the hcp- 
to-fcc hopping rate is at least 75 times higher than the 
fcc-to-hcp hopping rate, providing a lower bound for the 
energy difference of the two sites of | A£^| > kT ■ ln(75) = 
4meV. Alternatively, the hopping from hep to fee sites, 
which sets in at a temperature that is AT 1.5K lower 
than the onset for fcc-to-hcp hopping, allows a rough 
estimate of the energy difference for the two sites as 
AE Eb ' AT/T. Taking the experimental uncer- 
tainties into account, this provides an upper bound of 
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FIG. 1: By means of an extended site-mapping technique 
the adsorption site of Cu/Cu(lll) was determined to be the 
fee site. Individual adatoms were repeatedly repositioned 
using lateral manipulation at a temperature of 7 and 11 K. 
One STM image of a series and the corresponding model are 
shown. In the topmost part of the figure, the positions of 
these individual adatoms (circles at 7K, crosses at 11 K) as 
well as of a dimer and a trimer (arrows) with respect to the 
row are shown. 



A^ < 8meV, subject to the assumption that the fee and 
hep diffusion prefactors differ by no more than a factor of 
ten (in fact our theory calculations yield identical pref- 
actors of i/o = 1 X lO^^s-i [23). 

This result is consistent with a previous DFT study 
25] and congruent with our DFT calculations that yield 
an fee preference of about 6meV. An identical energy 
difference of 6 meV arises in a corresponding DFT study 
of monomer adsorption on a prerelaxed frozen surface, or 
when we increase the /c-point sampling while lowering the 
smearing. A similar fee preference arises when, instead, 
we use a four-layer 5-by-4 unit cell. These tests sug- 
gest that the monomer calculations are converged with 
respect to the choice of DFT parameters. The small 
monomer-adsorption preference reflects subtle differences 
in the interplay between the adsorbate-induced density 
reponse and the exact location of subsurface atoms, as 
we shall discuss further in a forthcoming publication. 

In experiments, dimer s do not exhibit any translational 
diffusion up to a temperature of 21 K. Nonetheless they 
already diffuse locally at a temperature as low as 7K, as 
can be seen from the unstable appearance of the dimer 
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FIG. 2: (a) At T = 7K Cu dimers appear unstable in STM 
images, showing that local diffusion takes place during scan- 
ning, (b) Long-range interactions with a monomer close by 
can stabilize the dimer. (c) At a temperature of 5K, the 
dimer diffusion is slow enough to trace individual diffusion 
steps. The white arrows indicate the fast scan direction. The 
crystal orientation corresponds to the one in Fig.Q] 



in the STM image shown in Fig.ISJa). As predicted by 
theory |l3] and verified by the site determination demon- 
strated in Fig.n the local diffusion takes place in a cell 
centered around one atom in the topmost layer of the 
substrate surface. This cell of local diffusion consists of 
three hep and three fee sites around the on-top position 
of the substrate atom (see scheme in Fig.^. In the fol- 
lowing we will use the notation //, hh, fh, fb, hb, and bb 
for fcc-fcc, hcp-hcp, fcc-hcp, fcc-bridge, hep-bridge, and 
bridge-bridge dimer occupation, respectively. 

As can be seen in Fig.[21^b) a dimer can be stabi- 
lized by long-range interactions 16, 26, 27] with another 
monomer positioned close to the dimer. Moreover, re- 
ducing the temperature to below 7K, the dimer intra-cell 
diffusion becomes so slow that it is possible to trace the 
individual diffusion steps by recording a series of STM 
images. The stable // and the metastable fh configura- 
tion can be observed (see Fig.I2fc) ). The relative prob- 
ability of observation at 5K {Pfh/Pff — 1/20) yields 
an energy difference of only AE = —kT ■ \ii{Pfh/Pff) = 
(1.3±0.5) meV. The hh configuration could not be found. 

An Arrhenius plot of the dimer // to fh hopping 
rate (Fiffjpidocuments that thermally assisted tunnel- 
ing [22, 12M enables the hopping of the dimers at very 
low temperatures (T < 5K), despite the large mass mcu 
of a Cu atom. In a simple model with a parabolic barrier 
at the transition state fb, characterized by an imaginary 
frequency a;/^ through E = Eb + 1/2 • the 
transmission probability due to tunneling for an energy 
e is known to be approximately exp(— 2i7r(£^B — ^)/^^fb) 
|30| . In the case of thermally assisted tunneling this prob- 
ability has to be multiplied by the probability for thermal 
activation to the energy e of exp(— e//cT). This results in 
a critical temperature Tc = hu;fi)/2i7Tk, which marks the 
limit between thermally activated hopping and thermally 
assisted tunneling. Above this temperature we obtain the 
conventional Arrhenius behavior with the activation en- 
ergy Eb^ whereas for the low-temperature regime, only 
thermally assisted tunneling to the ground state is rele- 
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FIG. 3: Arrhenius plot of the dimer // to fh hopping rate. 
The hopping rate is determined by direct thermal activation 
(process A) at temperatures above 6 K and by thermally as- 
sisted tunneling to the ground state (process B) below 5K. 
The dashed line represent a fit, taking only these two pro- 
cesses into account. At 4.5 K< T < 6K an even better fit 
(dotted line) is obtained by including also thermally assisted 
tunneling for all intermediate energies (processes C). Tunnel- 
ing at [/ = 10 mV, 7=10 pA excluded tip-induced diffusion. 
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FIG. 4: Illustration of the DFT analysis of the dimer intra- 
cell dynamics 1^ observed on Cu(lll). Solid (dashed) lines 
illustrate the variation in the fully relaxed dimer-configuration 
energies calculated using the 4 x 4 (2 x 2) /c-point sampling. 



vant, and the slope in the Arrhenius plot is determined by 
the energy difference AE of the two ground states. Only 
in the temperature region just below Tc do all intermedi- 
ate processes come into play (see fits in Fig.|3|. The fit to 
the data yields z/o = 8- 10^^^°-^ s"^ and Eb = 18±3meV. 
Assuming the barrier to be cosine-shaped, cofb and thus 
Tc can be estimated from the distance d = 1.47 A be- 
tween the f ee and the hep sites and the barrier height Eb 
as Tc — \/ Eb/^tticu ' h/kd = 6.0 K, in agreement with 
experiment. The surprising observation of the tunneling 
of a heavy Cu atom can be understood when considering 
that the atom not only has a very low barrier but also an 
exceptionally small distance to overcome. In contrast to 
the recently reported tunneling of carbon adatoms |2^ . 
we do not see any influence of descrete vibrational levels 
of the dimer. This indicates a very strong vibronic cou- 
pling of the dimer to the substrate, clearly ruling out the 
possibility of a Cu2/Cu(lll) quantum rotor 13]. 

From calculations using 2x2 /c-point sampling (Fig.^ 
dashed lines), we predict an intra-cell dimer-energy varia- 
tion that is in qualitative agreement with the original Al- 
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TABLE I: Site preferenees, kinetie barriers, and stability of 
monomers, dimers, and trimers on Cu(lll) as observed in 
low-temperature STM and ealeulated in first-prineiples DPT. 





O J. iV± 


X-ZJ? J. 


Monomer: 






pi Citli 1 CLI olut; 


fee 


fee 


diffusion barrier Eb (meV) 
attempt frequeney (s~^) 


37 ±5 


50 


5 X 10^^^^ 


1 X 10^2 


= ^hcp - Eicc (meV) 


4 < < 8 


6 


Dimer: 






preferred site 


fee-fee 


fee-fee 


intra- eell- diffusion Eb (meV) 


18 ±3 


28 


^fcc-hcp — ^fcc-fcc (meV) 


1.3 ±0.5 


21 


Row of adatoms: 






preferred site 


fee 


fee 



dimer study . As in these previous theoretical calcu- 
lations and consistent with the experiments, we find that 
intra-cell rotation remains possible even when low tem- 
peratures inhibit the net inter-cell diffusion, for which we 
determine a concerted-sliding barrier £^B,dif ^ 190 meV. 
Using the more accurate 4x4 /c-point sampling we pre- 
dict an energy variation (Fig.^ solid lines), in even bet- 
ter agreement with experimental observations. The 4x4 
/c-point calculations establish a near degeneracy between 
the //, //i, and hh dimer configurations and a signifi- 
cantly reduced effective rotation barrier, i. e. the energy 
of the fb configuration. The accuracy improvement ob- 
tained by using the 4x4 /c-point sampling in the dimer- 



configuration survey is selective and correlated with the 
dimer-induced surface distortion as also previously found 
in the Al-addimer study We attribute the experi- 

mental observation of dimer confinement to the rotation 
process, ff^fb^ fh (and back), which has a lower 
barrier than the f f ^ fb ^ fh ^ hb ^ hh rotation and 
the direct ff^bb^ hh transition. 

Trimers were found to be always of fcc-fcc-fcc type and 
stable in the temperature range of 5 — 21 K. 

In summary, we have shown that individual copper 
adatoms as well as adatoms in dimers, trimers, and larger 
structures prefer fee over hep sites on Cu(lll). The ex- 
perimentally observed difference in binding energy for the 
monomer is only about 6meV. Moreover, in the case of 
copper dimers, we were able to observe a confined intra- 
cell metal dimer motion for a wide temperature range 
of 4 — 21 K, which at low temperatures is enabled by 
thermally assisted tunneling. Even though the small en- 
ergy differences inherent in this system may attain the 
limits of current DFT accuracy, a detailed DFT analy- 
sis yields monomer and dimer energetics that are in very 
good agreement with the experiments (see Table P). 
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